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Solvent-swollen  polymer  gels  can  be  utilized  as  mechanical  simulants  of  biological  tissues  to  evaluate 
protective  systems  and  assess  injury  mechanisms.  However,  a  key  challenge  in  this  application  of 
synthetic  materials  is  mimicking  the  rate-dependent  mechanical  response  of  complex  biological  tissues. 
Here,  we  characterize  the  mechanical  behavior  of  tissue  simulant  gel  candidates  comprising  a  chemically 
crosslinked  polydimethylsiloxane  (PDMS)  network  loaded  with  a  non-reactive  PDMS  solvent,  and 
compare  this  response  with  that  of  tissue  from  murine  heart  and  liver  under  comparable  loading 
conditions.  We  first  survey  the  rheological  properties  of  a  library  of  tissue  simulant  candidates  to 
investigate  the  effects  of  solvent  loading  percentage,  reactive  functional  group  stoichiometry,  and  solvent 
molecular  weight.  We  then  quantify  the  impact  resistance,  energy  dissipation  capacities,  and  energy 
dissipation  rates  via  impact  indentation  for  the  tissue  simulant  candidates,  as  well  as  for  the  murine 
heart  and  liver.  We  demonstrate  that  by  tuning  these  variables  the  silicone  gels  can  be  engineered  to 
match  the  impact  response  of  biological  tissues.  These  experiments  inform  the  design  principles  required 
for  synthetic  polymer  gels  that  are  optimized  to  predict  the  response  of  specific  biological  tissues  to 
impact  loading,  providing  insight  for  further  tuning  of  this  gel  system  to  match  the  impact  response  of 
other  “soft  tissues”. 

©  2013  Published  by  Elsevier  Ltd. 


1.  Introduction 

The  defense,  law  enforcement,  automobile,  and  contact  sports 
industries  have  identified  a  growing  need  for  head  and  torso 
forms,  and  general  testing  media,  that  accurately  mimic  the  rate- 
dependent  mechanical  response  of  the  human  body  for  assess¬ 
ment  of  new  protective  equipment  and  for  understanding  injury 
mechanisms  in  response  to  a  broad  range  of  impact,  blunt  trauma, 
and  penetrating  threats  (Bertrand  et  al.,  2008;  Bresson  and  Frank, 
2010;  Roberts  et  al.,  2007;  Thai!  et  al.,  2002).  Common  simulants 
exploited  in  these  forms  for  compliant  tissues  comprising  skin, 
brain,  muscle,  and  internal  organs  include  ballistic  gelatin,  various 
rubbers,  leather,  silicone  elastomers,  soap,  lard,  and  clay  (Appleby- 
Thomas  et  al.,  2011;  Jussila  et  al.,  2005;  Merkle  et  al.,  2008). 
In  most  cases,  the  tissue  simulants  are  a  crude  representation 
of  the  biological  tissue  response;  however,  detailed  assessment  of 
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protective  equipment  and  injury  mechanisms  requires  more  accu¬ 
rate  simulants  that  quantitatively  mimic  complex  rate-dependent 
behavior  of  biological  tissues.  Therefore,  a  critical  challenge 
for  these  applications  is  to  develop  environmentally  stable, 
easily  processable,  and  cost  effective  synthetic  materials  that 
recapitulate  the  mechanical  response  of  biological  tissues 
(juliano  et  al.,  2006;  Moy  et  al.,  2006)  that  are  mechanically 
compliant  and  exhibit  complex  nonlinear,  time-dependent  defor¬ 
mation  behavior  (Bisplinghoff  et  al.,  2009;  Nicolle  et  al,  2010; 
Oliver  et  al.,  2010;  Saraf  et  al.,  2007;  Song  et  al.,  2007;  Storm  et  al., 
2005;  Weiss  et  al.,  2002). 

Many  different  gel  chemistries  offer  novel  methods  to  tune  gel 
properties  and  better  mimic  the  mechanical  response  of  biological 
tissues.  Examples  include  polybutadiene  (Lenhart  and  Cole,  2006), 
epoxy  (Mrozek  et  al.,  2012),  polyvinyl-alcohol  (Stammen  et  al., 
2001 ),  physically  associating  gels  (Seitz  et  al.,  2009),  and  various 
hydrogels  (Gong,  2010).  Block  copolymer-based  gels  containing  a 
mid-block  selective  solvent  have  been  explored  explicitly  for  use 
as  tissue  simulants  in  head  and  torso  forms  Quliano  et  al,  2006; 
Kalcioglu  et  al.,  2011 ;  Moy  et  al.,  2006).  The  mechanical  response 
of  these  materials  can  be  tailored  through  solvent  loading,  the 
incorporation  of  diblock  polymer,  the  ratio  of  the  constituent 
blocks,  and  the  manipulation  of  the  block  copolymer  chemistry. 
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However,  those  gels  synthesized  to  date  have  not  accurately 
mimicked  the  response  of  soft  tissues  such  as  heart  and  liver 
under  concentrated  impact  loading  (Kalcioglu  et  al.,  2011). 

Recently  developed  environmentally  durable  silicone  gels, 
composed  of  a  chemically  crosslinked  polydimethylsiloxane 
(PDMS)  network  loaded  with  a  non-reactive  PDMS  solvent,  offer 
the  potential  to  mimic  the  complex  mechanical  response  of  soft 
biological  tissues  (Mrozek  et  al.,  2011).  By  altering  the  concentra¬ 
tion  of  elastically  active  crosslinking  chains,  trapped  and 
untrapped  physical  chain  entanglements,  and  dangling  chain  ends, 
as  well  as  the  corresponding  molecular  weight  (MW)  of  these 
chains,  both  the  elastic  and  viscoelastic  properties  of  the  silicone 
gels  can  be  tuned  systematically  (Bibbo  and  Valles,  1984; 
Kawamura  et  al.,  2002;  Patel  et  al.,  1992;  Roth  et  al.,  2004; 
Urayama,  2008;  Vega  et  al.,  2001;  Villar  and  Valles,  1996).  Thus, 
in  addition  to  environmental  stability,  a  significant  advantage  of 
silicone  gels  is  the  ability  to  tailor  the  properties  by  altering  the 
network  structure,  solvent  type,  solvent  loading,  and  solvent  MW 
(Mrozek  et  al.,  2011 ).  Specifically,  when  the  solvent  MW  is  larger 
than  the  entanglement  MW  of  the  polymer,  such  PDMS-based  gels 
can  exhibit  significant  strain-rate  stiffening,  suggesting  the  poten¬ 
tial  to  incorporate  and  tailor  the  strain  rate-dependent  response 
for  tissue  simulant  applications.  However,  for  this  approach  to  be 
useful  for  the  design  of  tissue  simulants,  such  gels  must  be 
synthesized  to  match  the  approximate  stiffness  of  most  biological 
soft  tissues,  and  the  relevant  mechanical  response  of  interest— the 
extent  and  rate  of  energy  dissipation  under  impact  loading— must 
be  compared  directly  between  these  organogels  and  the  various 
hydrated  biological  tissues. 

To  that  end,  here  we  synthesized  and  analyzed  the  mechanical 
response  of  PDMS  gels  in  comparison  to  tissues  obtained  from 
murine  heart  and  liver.  Due  to  the  complexity  of  mechanical 
characterization  at  higher  rates,  rheological  experiments  at  low 
strains  and  low  shear  rates  were  first  conducted  to  provide  insight 
into  the  effects  of  solvent  loading,  solvent  MW,  and  network 
stoichiometry  on  the  frequency-dependent  mechanical  response 
as  compared  to  the  heart  and  liver  tissues.  Since  such  experiments 
cannot  quantitatively  predict  the  impact  response  of  the  materials, 
due  to  differences  in  loading  rates  and  configurations  as  well  as 
contact  geometries  (bulk  sample  response  vs.  local  response  in 
concentrated  impact  indentation),  next  we  performed  impact 
indentation  experiments  to  quantify  the  impact  resistance  and 
the  energy  dissipation  capacities  and  rates  of  these  gels  and 
tissues.  These  findings  provide  design  principles  and  approaches 
to  modulate  the  mechanical  response  of  tissue  simulant  gels  to 
concentrated  impact  loading,  as  required  for  the  development  of 
biofidelic  tissue  simulants. 


2.  Experimental 

2.1.  Materials 

Gel  processing  details  are  discussed  in  previous  work  (Mrozek  et  al.,  2011). 
These  gels  comprised  a  vinyl-terminated  PDMS  (v-PDMS;  reactive  network  pre¬ 
cursor)  and  methyl-terminated  PDMS  (m-PDMS;  non-reactive  theta  solvent).  Chain 
length  of  v-PDMS  was  117  kg/mol,  exceeding  the  entanglement  molecular  weight 
MWent  of  PDMS  (~29  kg/mol)  (Gent  et  al.,  1994)  and  leading  to  an  entanglement- 
dominated  network.  We  investigated  gel  properties  by  varying  the  solvent  loading, 
solvent  MW,  and  reaction  stoichiometry  (the  molar  ratio  of  silane  to  vinyl 
functional  groups).  The  stoichiometry  was  varied  from  4:1  (previously  determined 
to  provide  an  optimum  network  structure  (Mrozek  et  al.,  2011))  to  2:1.  Solvent 
loading  was  varied  from  10  vol%  to  80  vol%,  and  the  solvent  MWs  considered  were 
1.1, 139,  and  308  kg/mol. 

Heart  and  liver  organs  were  harvested  from  adult  rats.  All  experiments 
involving  animals  followed  the  University  lACUC  protocol  and  the  NIH  guidelines 
for  animal  care.  Tissue  discs  of  8  mm-diameter  and  thickness  of  3-5  mm  were 
prepared  using  a  surgical  punch,  and  all  tissues  were  stored  in  Krebs-Hensleit 
buffer  immediately  after  excision  and  throughout  all  experiments  reported  herein. 


2.2.  Rheology 

To  measure  the  shear  storage  modulus  C',  loss  modulus  G".  and  loss  tangent 
tan  S  (i.e.,  G"/G'  employed  as  a  measure  of  the  dissipation  of  deformation  energy 
(Ferry,  1980)),  rheological  measurements  on  PDMS  gels,  heart  and  liver  tissues 
were  conducted.  For  details,  see  Supplementary  Information  SI. 

2.3.  Instrumented  impact  indentation 

Impact  indentation  was  conducted  on  all  tissues  and  gels  with  a  stainless  steel 
flat  punch  probe  of  R=1  mm  at  25  “C  using  a  commercially  available  pendulum- 
based  instrumented  nanoindenter  (Micro  Materials  Ltd.,  UK).  The  development  of 
this  method  is  discussed  in  detail  in  previous  work  (Constantinides  et  al.  (2009), 
Constantinides  et  al.  (2008),  Kalcioglu  et  al.  (2011 )).  Impact  resistance  described  via 
maximum  penetration  depths  x„ax.  energy  dissipation  capacity  K,  and  energy 
dissipation  rate  described  via  a  quality  factor  Q  were  computed  from  the  acquired 
displacement  vs.  time  response  of  the  pendulum.  See  Supplementary  information 
S2  for  full  discussion  of  this  method  and  data  analysis.  Experiments  herein  were 
conducted  at  impact  velocities  ranging  0.4-2  cm/s,  corresponding  to  a  strain 
energy  density  of  2-20  kj/m^  comparable  to  other  macroscale  impact  testing 
methods  (Snedker  et  al.,  2005)  due  to  the  relatively  small  probe  contact  area  and 
volume. 


3.  Results  and  discussion 

3.1.  Comparison  of  gel  and  tissue  rheology  at  low  strains  and  rates 

3.1.1.  Effect  of  solvent  loading  on  the  magnitude  and  rate 
dependence  of  G',  G",  and  tan  S 

Solvent  has  a  significant  impact  on  the  modulus  of  these  gels  in 
two  ways;  (1)  the  solvent  will  dilute  the  polymer  network  to 
decrease  the  density  of  network  chains  per  unit  volume;  and 
(2)  the  solvent  will  dilute  the  number  of  trapped  entanglements. 
Both  of  these  effects  will  be  enhanced  with  increased  solvent 
loading.  Previous  work  with  these  silicone  gels  has  shown  that,  at 
50%  solvent  loading  the  low  frequency  G',  is  ~40  kPa  (Mrozek 
et  al.,  2011)  which  is  stiffer  than  many  biological  soft  tissues. 
Therefore,  silicone  gels  of  lower  stiffness  (greater  mechanical 
compliance)  are  needed  for  tissue  simulants  that  mimic  soft  tissue 
response.  Thus,  gels  were  synthesized  to  contain  10,  30,  40,  50,  60, 
70,  75,  80,  and  85  vol%  solvent  using  the  optimum  stoichiometry 
of  4;1  with  a  solvent  MW  of  1.1  kg/mol.  Note  that  the  optimum 
stoichiometry  was  determined  by  the  reaction  conditions  that 
provided  a  maximum  G'.  The  materials  cured  uniformly  at  solvent 
loadings  up  to  80  vol%;  at  85  vol%,  the  gel  consistency  was  not 
uniform  and  included  significant  volumes  of  uncured  precursor.  G' 
decreased  as  the  solvent  loading  increased  from  10%  to  80%  by 
approximately  two  orders  of  magnitude.  Further,  the  lack  of  strong 
frequency  dependence  in  this  measurement  range  was  consistent 
with  gels  containing  non-entangled  solvents  (Fig.  la).  When  G' 
measured  at  1  Hz  is  graphed  as  a  function  of  polymer  vol%  on  a 
log-log  scale,  the  data  are  well  fit  by  a  straight  line,  indicating 
power  law  behavior  (Fig.  Id).  The  obtained  scaling  factor  (slope  of 
the  linear  fit)  of  2.5  is  slightly  larger  than  the  theoretical  value  of 
2.3  for  a  network  formed  in  a  theta  solvent  (Obukhov  et  al.,  1994). 
G"  also  decreased  by  nearly  two  orders  of  magnitude  when  the 
solvent  loading  was  increased  from  10  to  80  vol%  (Fig.  lb).  This 
resulted  in  a  very  low  tan  S  in  the  range  of  0.07-0.005,  consistent 
with  low  capacity  for  energy  dissipation  (Fig.  Ic).  To  provide  a 
comparison,  small-strain  rheology  was  performed  on  murine  heart 
and  liver  tissue  (Supplementary  Fig.  S2).  The  approximate 
mechanical  response  of  these  tissues  is  denoted  in  Fig.  la-c,  and 
agree  reasonably  well  with  previous  macroscale  rheological 
experiments  for  liver  (G'~0.4-10  kPa  and  G"~0.1-0.4  kPa)  from 
various  animal  sources  (Georges  et  al.,  2007;  Kiss  et  al.,  2004;  Liu 
and  Bilston,  2000;  Ozcan  et  al.,  2011).  These  results  show  that 
although  both  the  magnitude  and  frequency  dependence  of  G'  of 
heart  were  reasonably  matched,  we  were  unable  to  match  the 
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Fig.  1.  Effect  of  1.1  kg/mol  solvent  loading:  (a)  storage  modulus  G',  (b)  loss  modulus  C",  and  (c)  loss  tangent  tan  S  with  silane:vinyl  stoichiometry  of  4:1  as  a  function  of 
frequency,  compared  to  tissue  from  murine  heart  and  liver,  (d)  G'  as  a  function  of  polymer  content  for  PDMS  gels  containing  1.1  kg/mol  solvent. 


frequency  dependence  of  G"  by  tuning  the  solvent  loading  %. 
Nevertheless,  this  comparison  indicates  that  these  gels  can  be 
tuned  to  exhibit  macroscale  storage  moduli  at  low  strains  and 
strain  rates  that  approximate  those  of  biological  soft  tissues. 


3.1.2.  Effect  of  solvent  MW  on  the  magnitude  and  rate  dependence 
of  C,  G",  and  tan  S 

Effects  of  solvent  MW  on  the  viscoelastic  properties  of  the 
present  gel  system  were  investigated  in  detail  previously  (Mrozek 
et  al.,  2011 ).  It  was  found  that  by  manipulating  the  solvent  MW  to 
effectively  modulate  the  number  of  entanglements  per  solvent 
chain,  the  gel  mechanical  response  could  be  tuned.  Gels  were 
synthesized  to  comprise  10,  30,  40,  50,  60,  70,  75,  and  80  vol% 
solvent,  using  the  optimum  stoichiometry  4:1  and  a  solvent  MW  of 
308  kg/mol.  Increasing  the  solvent  MW  from  1.1  kg/mol  to  308  kg/ 
mol  did  not  significantly  affect  G'  at  low  frequencies,  indicating  a 
consistent  network  microstructure  (Fig.  2a).  However,  higher 
solvent  MW  did  induce  an  enhanced  frequency  dependence  of 
G'  (Fig.  2a).  The  enhanced  frequency-dependent  G'  is  attributed  to 
the  longer  relaxation  times  of  the  solvent  chains  in  the  polymer 
network  as  the  solvent  MW  increases  (Doi  and  Edwards,  1986; 
Gennes,  1979).  Solvent  chains  are  expected  to  enhance  gel  stiffness 
when  the  measurement  time  scale  is  shorter  than  the  relaxation 
time  of  the  solvent  chains  (Urayama  et  al.,  2001),  producing  an 
increase  in  G'  at  higher  frequencies.  G"  also  increased  significantly 
for  solvent  of  higher  MW,  particularly  at  high  solvent  loadings 
(Fig.  2b).  This  is  attributed  to  the  ability  for  the  solvent  to  entangle 
with  the  network  and  with  itself,  providing  mechanisms  for 
energy  dissipation.  The  enhanced  energy  dissipation  results  in 
much  larger  tan  S  values  at  loadings  above  50  vol%  (Fig.  2c).  As  a 
result,  gels  containing  high  MW  solvent  are  too  energy  dissipative 


(larger  tan  5)  when  compared  to  murine  heart  and  liver  tissue  for 
gels  that  also  exhibit  the  relevant  range  of  storage  moduli  (i.e., 
solvent  loadings  of  70-80  vol%).  To  demonstrate  that  the  solvent 
MW  can  be  used  to  tailor  the  energy  dissipation  along  with  the 
strain-rate  dependence,  a  gel  was  synthesized  consisting  of  75  vol 
%  solvent  using  the  optimum  stoichiometry  ratio  and  a  solvent 
MW  of  139  kg/mol  (Fig.  2d-f):  this  gel  indeed  exhibited  a  low- 
frequency  G'  similar  to  that  of  the  gels  containing  1.1  and  308  kg/ 
mol  (Fig.  2d)  and  G"  and  tan  S  that  were  intermediate  to  the  higher 
and  lower  solvent  MWs  (Fig.  2e  and  f). 


3.3.3.  Effect  of  stoichiometry  on  magnitude  and  rate  dependence 
of  G',  G",  and  tan  S 

The  storage  modulus  of  the  gels  can  also  be  reduced  through  an 
increased  incorporation  of  defects  into  the  polymer  network 
structure.  Here,  we  examined  gels  at  a  range  of  crosslinker  ratios 
(4:1  to  2.25:1)  containing  60  vol%  solvent  with  a  MW  of  either 
1.1  kg/mol  or  308  kg/mol.  Fig.  2a-c  shows  that  decreasing  the 
stoichiometry  (for  the  1.1  kg/mol  MW  solvent  gel)  from  4:1  to 
2.25:1  reduced  G',  particularly  for  silane  to  vinyl  ratios  less  than 
3:1.  This  correlated  with  decreased  G'  at  low  frequencies  by  nearly 
an  order  of  magnitude,  which  was  similar  to  the  change  in  G' 
obtained  by  increasing  the  solvent  loading  from  60  to  80  vol%  at  a 
4:1  stoichiometry  (Fig.  la).  Despite  the  decrease  in  G',  G"  did  not 
change  significantly  as  a  function  of  stoichiometry  (Fig.  3b)  leading 
to  higher  tan  S  values  for  lower  crosslinker  ratios  (Fig.  3c).  Similar 
to  effects  of  solvent  MW  indicated  in  Fig.  2a-c  over  a  range  of 
solvent  loadings,  here  the  solvent  of  higher  MW  resulted  in 
enhanced  frequency  of  G'  and  in  increased  G"  and  tan  <5  over  a 
wide  range  of  stoichiometric  ratios  (Supplementary  Fig.  S3). 
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Fig.  2.  Effect  of  308  kg/mol  solvent  loading:  (a)  storage  modulus  C',  (b)  loss  modulus  G",  and  (c)  loss  tangent  tan  5  values  with  silane:vinyl  stoichiometry  of  4:1  as  a  function 
of  frequency,  compared  to  heart  and  liver.  Comparison  of  the  (d)  storage  modulus  G',  (e)  loss  modulus  G",  and  (f)  loss  tangent  as  a  function  of  frequency  for  PDMS  gels 
produced  with  a  silane:vinyl  stoichiometry  of  4:1  containing  75  vol%  solvent  of  1.1, 139,  and  308  kg/mol,  resepectively. 
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In  summary,  this  rheological  survey  of  a  large  library  of  PDMS 
gels  indicates  that  the  magnitude  of  G'  can  be  tuned  to  match  that 
of  soft  tissues  at  the  frequencies  studied  herein,  by  decreasing  the 
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Fig.  4.  Effect  of  solvent  loading  percentage  for  PDMS  gels  with  silane:vinyl  stoichio¬ 
metric  ratio  of  4:1  and  solvent  molecular  weight  of  1.1  kg/mol.  (a)  Maximum 
penetration  depth  x„ax  as  a  function  of  impact  velocity:  as  the  solvent  loading 
percentage  increases,  the  gels  become  less  impact  resistant:  this  effect  is  more 
pronounced  at  higher  impact  velocities,  (b)  Energy  dissipation  capacity  K  as  a  function 
of  impact  velocity:  K  increases  as  solvent  loading  percentage  increases:  this  effect  is 
more  pronounced  at  higher  impact  velocities,  (c)  Quality  factor  Q  as  a  function  of 
impact  velocity:  increasing  the  solvent  loading  percentage  decreases  Q  (i.e.,  increases 
the  energy  dissipahon  rate).  Data  are  represented  as  mean  +  standard  error,  and 
standard  error  bars  may  appear  smaller  than  the  symbols. 


silane:vinyl  stoichiometry  and  increasing  the  solvent  loading. 
When  the  solvent  MW  was  lower  than  the  entanglement  molecular 
weight  (MWent)  of  PDMS,  changing  the  solvent  loading  or  the 
stoichiometry  did  not  alter  the  frequency  dependence.  By  increas¬ 
ing  the  solvent  MW  to  308  kg/mol,  gels  were  obtained  that 
exhibited  an  increased  frequency  dependence.  While  matching 
the  approximate  stiffness  and  rate-dependent  small  strain  response 
of  biological  tissues  are  important  components  of  tissue  simulant 
development,  the  impact  resistance  and  energy  dissipation 
processes  at  higher  rates  and  larger  strains  are  critical  to  recapi¬ 
tulate  the  response  of  soft  tissues  to  abrupt,  concentrated  loading 
that  is  typical  of  intentional  or  accidental  impact  events. 

3.2.  Impact  indentation  experiments 

3.2.1.  Effect  of  solvent  loading  on  x„,ax,  K  and  Q 

Impact  indentation  experiments  were  conducted  on  heart  and 
liver  tissues  and  on  PDMS  gels  with  silane:vinyl  stoichiometry  of 
4:1,  solvent  MW  of  1.1  kg/mol,  and  a  range  of  solvent  loading 
including  60,  70,  75  and  80  vol%.  As  expected,  x^ax  increased  as 
the  solvent  loading  increased  (Fig.  4a),  indicating  that  the  gels 
become  less  impact  resistant  with  increasing  solvent  vol%.  This  is 
in  agreement  with  small-strain  rheology  (Fig.  la),  which  showed 
decreased  G'  with  increase  in  solvent  loading.  Representative 
stress  vs.  strain  responses  calculated  at  an  impact  velocity  of 
~0.7  cm/s  also  showed  that  the  maximum  forces  attained 
decreased  with  increased  solvent  loading  (Supplementary  Fig. 
S4a).  Increased  solvent  loading  also  led  to  an  increase  in  the 
impact  velocity-dependence  of  penetration  depths  (Fig.  4a)  and  in 
the  energy  dissipation  capacity  K  (Fig.  4b).  Finally,  energy  dissipa¬ 
tion  was  faster  (lower  Q)  for  higher  solvent  loadings  (Fig.  4c). 
Impact  resistance  of  heart  tissue  was  matched  most  closely  by  the 
gel  comprising  70  vol%  solvent.  However,  this  gel  dissipated  less 
energy  in  the  first  impact  event  (lower  K)  and  dissipated  this 
energy  more  slowly  (higher  Q)  than  both  heart  and  liver  tissues 
under  the  impact  velocities  considered.  Although  the  gel  compris¬ 
ing  80  vol%  solvent  most  closely  approached  the  impact  resistance, 
K  and  Q  of  liver  tissue,  all  of  these  gels  were  in  fact  more  impact 
resistant  and  exhibited  lower  capacities  and  rates  of  energy 
dissipation  than  liver  tissue.  Thus,  although  trends  as  a  function 
of  solvent  loading  are  clear,  modulation  of  only  this  parameter  is 
insufficient  to  mimic  the  impact  behavior  of  either  of  these  soft 
tissues. 


3.2.2.  Effect  of  solvent  MW  on  x^^xi  K  tind  Q_ 

Impact  indentation  experiments  were  conducted  on  PDMS  gels 
with  solvent  loading  of  80  vol%  and  silane:vinyl  stoichiometric 
ratio  of  4:1  for  solvent  MWs  of  1.1,  139  and  308  kg/mol.  Fig.  5a 
shows  that  impact  resistance  of  these  gels  increased  with  increas¬ 
ing  MW  of  the  solvent  (lower  x^ax)-  Further,  the  dependence  of 
Xmax  on  impact  velocity  was  enhanced  with  decreasing  solvent 
MW.  This  suggests  that  the  rate  of  stiffening  of  the  gels  decreased 
as  the  solvent  MW  decreased,  recalling  the  macroscale  rheological 
trends  of  decreased  frequency  dependence  of  storage  moduli  with 
decreasing  solvent  MW.  Such  increased  stiffness  and  impact 
resistance  are  consistent  with  the  concept  of  an  increased  number 
of  physical  entanglements  as  the  solvent  MW  surpasses  MWent- 
Energy  dissipation  capacity  (Fig.  5b)  and  rate  (Fig.  5c)  also 
increased  substantially  as  the  solvent  MW  surpassed  MWent. 
which  could  be  attributed  to  increased  physical  entanglements 
and  to  correspondingly  decreased  solvent  mobility.  In  fact,  gels  of 
solvent  MWs  of  139  kg/mol  and  308  kg/mol  dissipated  most  or  all 
of  the  energy  during  the  first  impact  leading  to  K~'l  and  Q  less 
than  unity  (dashed  line).  We  note  that,  although  adhesion  was 
mitigated  to  the  extent  possible  in  these  experiments,  gel 
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Fig.  5.  Effect  of  solvent  molecular  weight  on  impact  response  of  PDMS  gels  with  a 
4:1  precursor  chain  stoichiometry  and  80vol%  solvent  loading,  (a)  Maximum 
penetration  depths  as  a  function  of  impact  velocity:  gels  become  more  impact 
resistant  as  the  molecular  weight  (MW)  of  the  solvent  increases;  this  effect  is  more 
pronounced  at  higher  impact  velocities,  (b)  Energy  dissipation  capacity  K  as  a 
function  of  solvent  MW:  gels  of  larger  solvent  MW  are  more  dissipative  at  a  given 
impact  velocity,  (c)  Quality  factor  Q  as  a  function  of  impact  velocity:  increasing 
solvent  MW  increases  the  energy  dissipation  rate  (lower  Q),  to  the  extent 
that  for  gels  with  MW  >  1.1  kg/mol,  energy  was  dissipated  during  the  first 
impact  cycle  leading  to  Q  values  below  the  dashed  line.  Data  are  represented  as 
mean  +  standard  error. 


adhesion  to  the  impact  probe  surface  was  quaiitatively  greater  at 
higher  soivent  MWs;  this  could  serve  as  an  additional  mechanism 
for  energy  dissipation.  We  conclude  from  these  results  that  gels 


comprising  solvents  of  lower  MW  are  favored  candidates  to 
approximate  the  impact  response  of  heart  and  liver  tissues. 

3.2.3.  Effect  of  stoichiometry  on  Xmax.  K  tmd  Q 

To  investigate  the  effect  of  precursor  chain  stoichiometry 
(representing  the  crosslink  density  within  the  present  range)  on 
impact  responses,  PDMS  gels  loaded  with  solvent  of  1.1  kg/mol 
MW  at  60  vol%  were  considered;  stoichiometric  ratios  ranged  from 
2:1  to  4:1.  Fig.  6a  shows  that  impact  resistance  of  these  gels 
increased  for  a  given  impact  velocity  as  the  crosslinking  density 
increased  (i.e.,  as  stoichiometry  varied  from  2:1  to  4:1).  These 
trends  agreed  well  with  macroscale  rheology,  for  which  increasing 
crosslinking  density  resulted  in  increased  C'  (Fig.  3a).  Representa¬ 
tive  stress  vs.  strain  responses  calculated  at  an  impact  velocity  of 
~0.7  cm/s  also  showed  that  maximum  attainable  forces  increased 
with  increased  crosslinking  density  (Supplementary  Fig.  S4b).  Gels 
of  highest  crosslinking  densities  (3:1  and  4:1  ratio)  exhibited 
similar  x^ax.  indicating  that  further  increases  in  crosslinker  den¬ 
sity  would  be  unlikely  to  further  decrease  x^ax  because  the 
mechanical  response  is  already  entanglement  dominated;  the 
stoichiometric  ratio  was  already  optimal  at  4:1  (Mrozek  et  al., 
2011).  Increasing  the  amount  of  crosslinker  resulted  in  lower 
energy  dissipation  capacities  (lower  K,  Fig.  6b)  and  lower  energy 
dissipation  rates  (higher  (J,  Fig.  6c),  due  ostensibly  to  an  increase  in 
the  elastic  chains  and  a  decrease  in  the  incorporation  of  dangling 
chain  ends.  Impact  resistance  and  energy  dissipation  rate  of  heart 
tissue  were  best  matched  by  gels  comprising  60  vol%  of  1.1  kg/mol 
solvent  at  a  precursor  stoichiometric  ratio  of  2.75:1.  Impact 
resistance  of  liver  tissue  was  best  matched  by  gels  comprising 
60  vol%  of  1.1  kg/mol  solvent  at  a  precursor  stoichiometric  ratio  of 
2:1,  though  this  candidate  tissue  simulant  dissipated  more  energy 
than  liver  tissue  (higher  K). 

3.2.4.  Comparison  of  impact  response  among  gels  and  soft  tissues 

The  above  discussion  makes  clear  that  the  three  design  para¬ 
meters  considered  (solvent  vol%,  solvent  MW,  and  precursor  ratio) 
can  result  in  coupled  or  competing  effects  on  the  impact  perfor¬ 
mance  metrics  of  interest  (x^ax.  K  and  Q).  Quantification  and 
understanding  of  this  coupling  allows  us  to  move  from  many 
iterative  experiments  to  reasonably  predictive  design  of  impact 
performance.  Here,  we  demonstrate  the  tuning  of  these  “knobs”  to 
realize  a  tissue  simulant  for  heart  tissue  impact  response,  and 
discuss  how  these  same  principles  can  be  extended  toward  design 
of  mechanical  tissue  simulants  for  other,  more  compliant  tissues 
such  as  liver. 

Heart:  Our  impact  indentation  experiments  showed  that  two 
gels  exhibited  similar  impact  resistance  (i.e.,  x^ax)  to  that  of  heart 
tissue  for  the  impact  velocities  studied  here  (Fig.  7a).  These  gels 
comprised  solvent  of  1.1  kg/mol  MW,  with  loading/precursor  chain 
stoichiometry  of  70vol%/4:l  or  60  vol%/2.75:l.  In  these  gel 
families,  we  observed  that  impact  resistance  could  be  tuned  by 
opposing  design  parameters:  increasing  the  solvent  loading  from 
60  vol%  to  70  vol%  increased  Xmax.  and  increasing  the  crosslinking 
density  (changing  the  stoichiometry  from  2.75:1  to  4:1)  decreased 
Xmax-  Although  these  two  gels  exhibited  similar  impact  resistance, 
the  candidate  with  the  lower  solvent  loading  and  crosslinking 
density  dissipated  more  energy  during  the  first  impact  event 
(larger  K),  and  dissipated  this  energy  more  quickly  (lower  Q) 
(Fig.  7b  and  c).  Thus,  this  tissue  simulant  candidate  better  matches 
the  impact  response  of  heart  tissue  in  terms  of  the  magnitudes  of 
Xmax.  K.  and  (2,  and  the  impact  velocity-dependence  of  Xmax  and  K. 
More  importantly,  by  tuning  the  crosslinker  density  and  solvent 
loading,  we  individually  tuned  the  three  metrics  of  impact 
behavior  to  better  mimic  the  response  of  a  specific  biological 
tissue.  Finally,  although  this  gel  exhibited  a  slight  decrease  in  the 
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Fig.  6.  Effect  of  silaneivinyl  stoichiometric  ratio  for  PDMS  gels  with  60  vol%  1.1  kg/ 
mol  molecular  weight  solvent,  (a)  Maximum  penetration  depths  as  a  function  of 
impact  velocity:  impact  resistance  of  gels  increases  (lower  Xmax)  by  increasing  the 
crosslink  density;  this  effect  is  more  pronounced  at  higher  impact  velocities. 
Impact  resistance  of  heart  and  liver  tissues  is  matched  by  the  PDMS  gels, 
(b)  Energy  dissipation  capacity  K  as  a  function  of  impact  velocity:  gels  dissipate 
more  energy  at  lower  crosslinking  densities.  Energy  dissipation  capacity  of  the 
tissues  is  matched  by  the  gels  with  lower  crosslinking,  (c)  Quality  factor  Q  as  a 
function  of  impact  velocity:  gels  dissipate  energy  more  quickly  (lower  Q)  as  the 
crosslinking  density  decreases.  Both  heart  and  liver  show  comparable  energy 
dissipation  rates  to  the  gels  with  lower  crosslinking  density.  Data  are  represented 
as  mean  ±  standard  error,  and  error  bars  may  appear  smaller  than  data  symbols. 


magnitude  of  Q  at  higher  impact  velocities,  Q.  of  heart  did  not 
depend  on  impact  velocity  for  the  velocities  explored  herein. 
Based  on  our  findings,  one  could  speculate  that  by  decreasing 
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Fig.  7.  (a)  By  changing  the  solvent  loading  vol%  and  silane:vinyl  stoichiometry,  we 
compare  two  different  gels  that  both  mimic  the  impact  resistance  of  heart  tissue  for 
the  impact  velocities  explored  herein,  (b)  These  two  gel  systems  exhibit  distinct 
energy  dissipation  capacities  K  (b)  and  quality  factor  Q(c).  For  the  impact  velocities 
considered,  the  gel  with  60  vol%  solvent  of  1.1  kg/mol  MW  and  silane:vinyl 
stoichiometry  of  2.75:1  most  closely  mimics  heart  tissue  by  all  three  metrics  of 
mechanical  impact  response. 


the  solvent  loading  vol%  to  below  60%,  the  velocity-dependence  of 
Q  in  these  gels  would  decrease  (Fig.  4c),  and  better  match  that  of 
heart  tissue.  However,  this  change  would  also  lead  to  an  increase 
in  the  magnitude  of  Q,  and  a  decrease  in  the  magnitude  and 
velocity-dependence  of  x^ax  and  K,  resulting  in  a  mismatch  of 
those  parameters.  To  compensate  for  this  outcome,  the  silane:vinyl 
stoichiometric  ratio  may  be  slightly  decreased  to  below  2.75:1, 
which  would  cause  an  increase  in  the  magnitude  of  Xmax  and  K, 
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Fig.  8.  (a)  The  two  top  candidates  for  matching  liver  tissue  impact  resistance  are 
the  PDMS  gels  with  60vol%  solvent/silane:vinyl  stoichiometry  of  2:1,  and  80vol% 
solvent/silaneivinyl  stoichiometry  of  4:1,  both  with  solvent  MW  of  1.1  kg/mol.  The 
two  top  candidates  for  matching  energy  dissipation  capacity  K  (b)  and  quality 
factor  Q  (c)  of  liver  tissue  are  the  PDMS  gels  with  60  vol%  solvent/silane:vinyl 
stoichiometry  of  2.25:1,  and  60vol%  solvent/silane:vinyl  stoichiometry  of  2.5:1, 
both  with  solvent  MW  of  1.1  kg/mol. 


a  decrease  in  the  magnitude  of  Q,  and  an  increase  in  the  velocity- 
dependence  of  Xmax  without  significantly  affecting  the  velocity- 
dependence  of  Q  (Fig.  6). 

Liver:  We  identified  several  gels  that  exhibited  either  x^ax 
(Fig.  8a),  K  (Fig.  8b),  or  Q  (Fig.  8c)  that  were  similar  to  those  of 
liver  tissue  for  the  impact  velocities  studied  here.  However,  unlike 
for  the  case  of  heart  tissue,  none  of  these  gels  captured  the 


magnitude  and  velocity-dependence  of  impact  response  of  liver 
tissue  by  all  metrics.  For  example,  gels  with  a  composition  of 
80  vol%—4:l-l.l  kg/mol  and  60  vol%—2:l-l.l  kg/mol  exhibited 
Xmax  similar  to  liver  tissue,  while  the  former  dissipated  less  energy 
and  the  latter  dissipated  more  energy  than  liver.  Further  increas¬ 
ing  the  MW  of  the  80%  4:1  gel  from  1.1  to  139  kg/mol  increased  the 
magnitude  of  K  but  also  decreased  the  magnitude  of  x^ax  and 
velocity-dependence  of  K  significantly.  On  the  other  hand,  gel 
systems  with  a  composition  of  60  vol%— 2. 5: 1-1.1  kg/mol  and 
60  vol%— 2.25:1-1.1  kg/mol  exhibited  similar  magnitude  of  K 
values  to  liver  but  yielded  lower  magnitude  of  x^ax  and  a  decrease 
in  the  impact  velocity-dependence  of  K.  Although  It  is  beyond  the 
scope  of  the  current  study  to  further  optimize  these  gels  toward 
liver  tissue  simulants,  we  can  draw  on  the  design  parameters 
identified  above  to  now  suggest  further  iterations  toward  that  aim. 
Based  on  the  gels  considered  herein,  an  improved  liver  tissue 
simulant  could  plausibly  be  attained  by  increasing  the  solvent  vol% 
in  gels  with  low  crosslinking  density  (between  60%  and  70%  for 
2.25:1  and  2.5:1  gels);  this  tuning  would  be  expected  to  Increase 
the  magnitude  of  x^ax  and  the  velocity-dependence  of  x^ax  and  K. 
However,  this  increase  in  the  solvent  vol%  would  also  likely 
increase  the  magnitude  of  K  values  to  be  higher  and  decrease 
the  magnitude  of  Q  values  to  be  lower  than  that  of  liver.  To 
compensate  for  this  change,  crosslinking  density  could  be 
increased  only  slightly  (precursor  ratios  between  2.5:1  and 
2.75:1)  to  achieve  lower  K  and  higher  Q  values  without  signifi¬ 
cantly  altering  the  magnitude  of  x^ax  and  the  velocity-dependence 
of  K,  0,  and  x^ax- 

In  summary,  we  Identified  design  principles  required  to  simul¬ 
taneously  optimize  key  metrics  of  impact  response  for  PDMS-based 
gels.  We  thus  demonstrated  specific  gel  compositions  that  matched 
the  impact  response  of  heart  tissue  with  reasonable  fidelity.  These 
findings  inform  the  design  of  synthetic  gels  to  recapitulate  the 
response  of  specific  biological  materials  under  impact. 
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